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Abstract 

The strange quark mass is calculated from QCD sum rules for the divergence of the vector 
as well as axial-vector current in the next-next-to-leading logarithmic approximation. The 
determination for the divergence of the axial-vector current is found to be unreliable due 
to large uncertainties in the hadronic parametrisation of the two-point function. 

From the sum rule for the divergence of the vector current, we obtain a value of m s = 
m s (l GeV) = 189 ± 32MeV, where the error is dominated by the unknown perturbative 
0{pc Z g) correction. Assuming a continued geometric growth of the perturbation series, we 
find m s = 178 ± 18MeV. Using both determinations of m s , together with quark-mass 
ratios from chiral perturbation theory, we also give estimates of the light quark masses m u 
and irtd- 
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1 Introduction 



Quark masses are amongst the fundamental set of "a priori" unknown parameters in the 
standard model (SM) of particle physics. They are not physical observables, as there do 
not exist free quarks in nature,^ and hence depend on the renormalization prescription 
applied to the model, but have the same status as the strong QCD coupling constant a s . 

In this work, we shall be concerned with a determination of the strange quark mass, 
m s , at the next-next-to-leading order (NNLO) in the framework of QCD sum rules [|1]. 
(For an overview see 0-fl-) Already a considerable number of determinations of the 
strange quark mass can be found in the literature |]5|-[T9|, however partly being incompatible 



within their errors. For this reason, and also because in the meantime there has been 
progress both on the experimental input, as well as on the theoretical expressions, we find 
it justified to reconsider the evaluation of m s . In addition, there is great interest in a 
precise value of m s for the calculation of direct CP-violation in the Kaon system within 
the SM, because the dominant matrix elements of four-quark operators contributing to 
e'/e scale like m' 2 



The basic object which is investigated in the simplest version of QCD sum rules is the 
two-point function ^(q 2 ) of two hadronic currents 

^(q 2 ) = i Jdxe^(n\T{j r (x)j r (0)^\n), (1.1) 

where Q denotes the physical vacuum. To be specific, in the rest of this work jr{x) will be 
the divergence of vector and axial- vector current, j(x) and js(x) respectively; 

j(x) = i(M — m) :Q(x)q(x): , j 5 (x) = i (M + m) : Q(x) 75 q(x) : , (1.2) 

with M and m being the masses of Q(x) and q(x). Note that these currents are renor- 
malization invariant operators. Throughout the paper, for notational simplicity, we shall 
drop the index for the pseudoscalar current, but we keep differing signs where they appear. 
The upper sign will always correspond to the divergence of the vector and the lower sign 
to the axial-vector current. It should be obvious that ^{q 2 ) is simply proportional to the 
two-point function of scalar (pseudoscalar) currents, 

^(g 2 ) = (MTm) 2 Tl s > p (q 2 ), (1.3) 



where we adopted the notation of ref. [22 



except presumably for the top quark 
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After taking two derivatives of ^{q 2 ) with respect to q 2 , ^"(q 2 ) vanishes for large q 2 , 
and satisfies a dispersion relation without subtractions (for the precise conditions see |23| ) : 



3C 



*^ = w? nq) = 2 hs-7-icY ds ' (L4) 





where p(s) is defined to be the spectral function corresponding to \P(s), 

p(s) = -Im*(s + ie). (1.5) 

7T 

To suppress contributions in the dispersion integral coming from higher excited states, 
it is further convenient to apply a Borel (inverse Laplace) transformation to eq. ( |1.4| ) which 
leads to PI 

B u V'(q 2 ) = = - f e- s / u p(s)ds, (1.6) 

u 6 J 


or slightly rewriting eq. ( p..6|) we obtain 

u^"(u) = J e~ s/u p(s)ds. (1.7) 


As we shall discuss in detail below, the left-hand side of this equation is calculable in 
renormalization group improved perturbation theory, if u can be chosen sufficiently large. 
Because the Borel transformation removes the subtraction constants in the dispersion re- 
lation and satisfies the identity 

B u ^(q 2 ) = u 2 B u V'(q 2 ), (1.8) 

we could have worked directly with ^(q 2 ). To be able to investigate other types of sum 
rules in the future, we nevertheless prefer to express our sum rule in terms of \l/"(g 2 ). 

Under the crucial assumption of quark-hadron duality, the right-hand side of eq. ( |1.7| ) 
can be evaluated in a hadron-based picture, still maintaining the equality, and thereby 
relating hadronic quantities like masses and decay widths to the fundamental SM parame- 
ters. Generally, however, from experiments the phenomeno logical spectral function p p h(s) 
is only known from threshold up to some energy s - Above this value, we shall use the 
perturbative expression pth{s) also for the right-hand side. This is legitimate if sq is large 
enough so that perturbation theory is applicable. The central equation of our sum-rule 
analysis then is: 

u 3 ^ h (u) = Je- s / u Pph (s)ds + Je- S / U p th (s)ds. (1.9) 

S 



3 All relevant formulae for the Borel transformation are collected in appendix A. 
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In addition, in our analysis we shall use the derivative of this equation with respect to u 
( "first-moment sum rule" ) : 

so oo 

= / - e- s/u p ph (s) ds + [- e- s ' u Pth (s) ds . (1.10) 

J Hi J %L 

so 

In sect. 2, we give the expressions for the theoretical part of the sum rules for scalar and 
pseudoscalar currents at the next-next-to- leading order, and in sect. 3, the phenomeno- 
logical parametrisations of the two-point functions are discussed. Sect. 4 deals with the 
perturbative continuum contribution and in sect. 5, we perform the numerical analysis of 
the sum rules. Finally, in sect. 6, we compare our results to previous determinations of m s 
published in the literature and give estimates for the light quark masses m u and m^. 



u 



a 

du 



u 



th 



U 



2 The theoretical two-point function 



In the framework of the operator product expansion (OPE) |24| the two-point function 
(1.1) can be expanded in inverse powers of Q 2 = —q 2 : 



*(Q 2 ) = (M T m) 2 g 2 {vI/ + v[; 2 i- + v[; 4 i- + v[; 6 i- + ...|. (2.1) 

The ty n contain operators of dimension n, and their remaining Q 2 -dependence is only log- 
arithmic. In this work, we shall treat mass like an operator of dimension 1. Otherwise, the 
contribution would be absent because one cannot construct a gauge- invariant operator 
of dimension 2. Likewise, ^"(Q 2 ) has the expansion 

*" (g2) = { *o + + K+i + K^S + ■■■}■ (2-2) 

In the following sections, we shall calculate these expressions explicitly. 

2.1 The perturbative contribution 

In general the perturbative expression for ^q(Q 2 ) is given by 



AT OO 1+1 

MQ 2 ) = ^E^E%^ (2-3) 

" i=0 j=0 



where iV is the number of colours, a = a s /ir and L = ln(Q 2 /[i 2 ). fi 2 is a renormaliza- 
tion scale. Throughout this work, we shall use the modified minimal subtraction scheme 
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MS |p5|| . From eq. ( |2.3| ), we easily obtain the corresponding expression for ^q(Q 2 



AT OO 1+1 

<(Q 2 ) = ^E^E<w 

i=0 j=l 



+ (j _ 



-2 



(2.4) 



Let us remark that the constant coefficients Cjo have dropped out of this expression. 

Because ^"(Q 2 ) is related to a physical quantity — the spectral function p(s) — it has 
to satisfy a homogeneous renormalization group equation (RGE): 



f d d 



(M T m) 2 ^(Q 2 ) = 0. 



(2-5) 



where (3(a) is the QCD /5-function and 7(a) is the mass anomalous dimension. In ap- 
pendix B, we have collected the coefficients of (3(a) and 7(a) in our notation. In eq. (|2.5| ), 
all coefficients multiplying a certain term a l U have to vanish separately. From this condi- 
tion one derives relations between the various Cjj which are also given in appendix B. The 
independent coefficients are conveniently chosen to be the Cn. The coefficients Cqi and en 
can be straightforwardly calculated and are found to be 



and 



Cll 



— C f . 
4 



(2.6) 



with Cf = (N 2 — 1)/2N. The NNLO coefficient C21 has been calculated recently in ref. |26 
The result is 



C21 



N=3 



10801 



(lr-i c(3) l^ 



893 31 



65 
24 



64 



rC(3) / 



— --C(3) )N- ---C(3) / 



65 



1 



32 



(2.7) 



where / is the number of flavours, and ((z) is the Riemann (^-function. Since VP" satisfies a 
homogeneous RGE, we can sum up the logarithms by choosing fi 2 = Q 2 . Then the coupling 
and the masses become running quantities, evaluated at Q 2 . To NNLO and N = f = 3 
this yields for \&q': 



N 



( 



5071 35 



yC(3))a(Q 2 ) 2 + 0(^ 



(2.* 



V 144 

The final step in the evaluation of the perturbative contribution to the two-point func- 
tion consists in performing the Borel transformation. Unfortunately, the running coupling 
a(Q 2 ) as well as the running masses M(Q 2 ) and m(Q 2 ) at NNLO are now complicated 
functions of Q 2 , which are conventionally expanded in powers of In -1 (Q 2 / A 2 ). The Borel 
transform of the resulting expressions cannot be given in closed form, but has to be ex- 
panded in powers of ln -1 (w/A 2 ), or calculated numerically P7|,p5 . 
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There is a different way to obtain the Borel transform of the two-point function. Because 
the differential operators for the RGE and the Borel transform act on different variables, 
/i 2 and Q 2 respectively, they commute PD[ . We thus apply the Borel operator directly to 
eq. ( |2.4| ). The necessary formulae are given in appendix A. One can then easily convince 
oneself that the resulting expression for (M =p m) 2 ^Q(u) again satisfies a homogeneous 
RGE, and the logarithms can be summed up through the choice fi 2 = u. This results in 
the running coupling and masses being evaluated at the scale u. Our final expression for 
§o(m) is 



K u) 



N 



l + a(tt)(cn+2(l-7 B )c 12 ) 



■ov>> - g7r2 

7T 2 

+ a 2 (c 21 + 2(1 - lE )c 22 + (3t| - 6 7B - y)c 23 ) (2.9) 

/ 7T 2 

+ a 3 (c 31 + 2(1 - lE )c 32 + (3 7 1 - 6 lE - — )c 33 - 2(2 7 | - 6 7 | - 7£ 7r 2 + vr 2 + 4C(3))c 34 



where 7 # is the Euler constant. Here, analogously to eqs. (|2.1|) and (|2.2|) , ty"(u) has been 
expanded in inverse powers of u: 

<?.,, . (M(w) =F m(u)) 2 f -„ -.,1 1 ~„ 1 1 . 

*"(«) = 1 1 J ^ 1 JJ { *S + + ^ + + • • • ) • (2-10) 

In our numerical analysis we have verified explicitly that our treatment of the Borel trans- 
form leads to the same result as the method of refs. |27,|28|, up to differences which are of 
higher order. 

In addition to eq. fl2.9|) , for the first-moment sum rule of eq. Q1.10 ) it is convenient to 
define the function 

= u^-$"{u). (2.11) 
ou 

Expanding this function equivalently to *$f"(u), eq. (|2.10|) , we obtain 



A f l + 5(«)(cii - 2 7b ci 2 ) + a 2 (c 21 - 2 7B c 22 + (3 7 | - 6 - y)c 23 ) (2.12) 



8tt 2 

+ a 3 (c 31 - 2 7B c 32 + (3 7 | - 6 - ^ )c 33 - 2(2 7 | - 12 lE - 1e ti 2 + 4C(3))c 34 ) 



_2 



To be able to estimate the uncertainty coming from the as yet unknown 0(a 3 ) term in 
the perturbative result, we have included this contribution which depends on the coefficient 
c 3 i in the general expressions of eqs. ( |2.9| ) and ( |2.12[ ), and in our numerical analysis we 
varied this coefficient in order to simulate a contribution from higher orders. Let us point 
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out however, that the 0(a 3 ) correction cannot be included in a completely consistent way, 
because this would also require the four-loop coefficients of the /3-function and the mass 
anomalous dimension, which have not yet been calculated. 



2.2 Dimension 2 operators 



The contribution to the dimension 2 operators can be obtained by expanding the pertur- 
bative result for the vector-current two-point function by Generalis |5D] in powers of the 
quark masses. Then the expression for the scalar (pseudoscalar) two-point function is cal- 
culable from the Ward-identity between the vector (axial- vector) and scalar (pseudoscalar) 
two-point functions (see e.g. eq. (2.3) of ref. f2"2]| ). There is, however, a subtlety, because 
this Ward-identity involves the renormalized, non-normal-ordered quark condensate (ipip). 
Since from a straightforward application of perturbation theory we get expressions con- 
taining the normal-ordered condensate (: ipip :), we still have to subtract from Generalis 
result a contribution which stems from the relation between (ipip) and (-.tpip:) |31j| . This 
additional contribution removes mass logarithms which are present at intermediate steps 



of the calculation p!6|,|22]]. For ty 2 we find 

* 2 (Q 2 ) 



N 



8tt 2 



± 



2L-2 + aC F (-3L 2 + 8L- y + 6C(3)) (ilf 2 + 
2L-4 + aCV(-3L 2 + 14L-22 + 6C(3)) Mm j 



(2.13) 



This result is in agreement with Chetyrkin et al. [32], who also performed a calculation for 
the vector current, but already gave expressions in terms of non-normal-ordered conden- 
sates. 



From eq. (|2.13j ) we obtain 
N 



4vr 2 



| l+aC F (-3L+7) (M 2 +m 2 )± l+aC F (-3L+lo) Mm J . (2.14) 



Proceeding for the Borel transform as in section 2.1, we finally find 

N 



and 
$ 2 (n) 



4tt 2 



N 
2^2 



| 1 + 0^(4 + 37 B ) (M 2 + m 2 )± l + aC F {j + ^ E ) Mm\, (2.15) 



i 
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l+aC F l^—+3j E 



(M 2 +m 2 )± 



Mm 



(2.16) 



Again, ^(u) satisfies a homogeneous RGE which justifies the exchange of Borel transfor- 
mation and renormalization group summation. 
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2.3 Dimension 4 operators 



There are three sources of dimension 4 operators contributing to ^(Q 2 ): explicit mass 
corrections ~ m 4 , terms originating from the quark condensate ~ m (?/>?/>) and the gluon 
condensate contribution ~ (aFF). To the order considered,^ the explicit quark mass as 
well as the gluon condensate contribution can, for example, be taken from ref. |p2 |. As 
for the quark condensate, to NNLO we also need 0(a s ) corrections. These can either be 
calculated directly |$3| , or inferred from the cancellation of mass logarithms after expressing 
('■^i^') through (ipip) [pT6|| , yielding identical results. The explicit expressions are: 



4,m 4 



N 
16vr 2 

1 



M 4 + 4M 2 m 2 + m 4 + 2(M 4 ± 2M 3 m ± 2Mm 3 + m 4 )L 



(2.17) 



2 
± 
1 



3 11 
1 + aC F [ --L + — 



4,FF 



1 + aC F --L + - 



(aFF) . 



(M(QQ)+m(qq)) 
(m(QQ) + M(qq) 



(2.18) 
(2.19) 



From these expressions one can immediately obtain ^>'l- ft * s an instructive exercise to 
show that again (M =p m) 2 ^' does satisfy the homogeneous RGE of eq. ( |2.5| ). To this 
end, the condensates are conveniently rewritten into renormalization group invariant con- 
densates |31~1, |34|1 . The formulae for the renormalization group invariant condensates are 
collected in appendix C. In addition, to show this one needs the terms of order a s m A which 
can also be calculated via expanding the result by Generalis |30| . 

As before, the next steps are performing the Borel transformation and summing up the 
logarithms through the choice fi 2 = u. This leads to 



+ 



^4 + ^aE^-^Em 4 



1 + aC F l — + -7 E 



Wi /ll 3 . 

— + Wi \C F I — + - J E ) + ^2 



± 



2 + aC F (7 + 3 lE ) (n™ Q + n A / q ) 



M +m 



± 



2^ + wi \C F (7 + 3 7£ ) + 2w 2 ) + 7iV< 



(M 3 m + Mm 3 



^M 2 m 2 



[2.20) 



For our counting of orders in perturbation theory see the remarks in appendix B. 
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and 



+ 



_!/lo4+?-av^r 9i -— 

2 \ 4 A Y 4 Y 

(flf Q + ± [2 + aC F (8 + 3 7 s) 1 (ft™ Q + ftf q ) 

(m 4 + 



, , 13 3 

1 + aC F ( — + -7e 



a 



— + wi[ C F [ — + -7s 



/13 3 



± 



2^ + wi (C F {8 + 3j E ) + 2w^j + 7° (- + 7s 



(M 3 m + Mm 3 ) ] 1(2.21 



with f2 3 and being the renormalization group invariant condensates as defined in ap- 
pendix C. All other definitions can be found in appendices B and C. 



2.4 Dimension 6 operators 

For the dimension 6 operators, we take into account only the most important contributions. 
In addition, since until today the renormalization group behaviour of dimension 6 operators 
has not been exploited completely, we shall neglect the running of those operators. Anyhow, 
numerically their contribution is only a very small correction, due to the suppression by 
powers of Q 2 . This will be discussed further in our numerical analysis of section 5. 
Including operators with up to one power in the quark masses, we have: 



± 



M(gqaFq) + m(gQaFQ) ± A^ 2 a{Qa^ v t a q qa^t a Q) 



4 



+ ^r: 2 a(Q lfl t a Q + q lf ,t a q ]T q A ^t a q A ) . 

A=u,d,s 



(2.22) 



Here, (gipaFip) is the so called "mixed" condensate, and the other two terms are four-quark 
operators. Since there exist no reliable estimates for the vacuum expectation values of these 
four-quark operators in the literature, we follow the usual procedure in sum rule analyses 
by using the vacuum dominance hypothesis [|I|] to relate them to the quark condensate. 
After taking two derivatives and performing the Borel transformation this leads to 



± 



1 



M(gqoFq)+m(gQoFQ) 



4C 



3iV 



F 2- 

7T a 



(QQ) 2 + (qq) 2 ±9(QQ)(qq) . (2.23) 



The corresponding quantity $6 for the first moment sum rule is simply given by — 4\?g. 



S 



3 Phenomenological two-point function 



In this section, we have to distinguish between the scalar and the pseudoscalar two-point 
function. Let us therefore discuss both separately. In addition, we now have to specify the 
actual flavour content of the currents of eq. ( |1.2| ). Since in the end we want to calculate 
the strange quark mass, Q should be chosen to be the strange quark and M = m s . The 
other quark is conveniently chosen to be a light quark, because for these channels we have 
reasonably good information on the spectral functions. For defmiteness, we shall use the 
up-quark and m = m u , neglecting isosp in-breaking effects. 

3.1 The scalar two-point function 

Generally, the phenomenological spectral function p(s) is given by 

pr{s) = (2n)^(n\j r (0)\T)(T\j r (0y\n)5\q-p r ), (3.1) 

where T are intermediate states with the correct quantum numbers over which we have to 
sum and calculate the corresponding phase-space integrals. 

For the scalar two-point function, the lowest lying state which contributes to the spectral 
function is the i^7r-system in an s-wave I = 1/2 state. The contribution of this intermediate 
state yields the inequality || 

3 1 

Pph (s) > 0(s - s+ ) — - yfc - s + )(s - S _) \d(s)\ 2 , (3.2) 

where 

s + = (M^+ + M w o) 2 , s_ = (M K + — M^o) 2 , (3.3) 
and d(s) is the strangeness changing scalar form factor which appears in Ki 3 decays: 

(AP')\ (Sl,u)(0)\K + (p)) = -±={(p + p%f + (s) + (p-pV-00] , (3.4) 

and 

d(s) = (M 2 K -M 2 n )f + (s) + sf4s). (3.5) 

The physical region for decays is mf < s < s_, whereas in the spectral function s > s + . 
In this region, d(s) is not directly accessible to experiment. 

In the following, we assume that the spectral function is saturated by a sum of Breit- 



Wigner resonances multiplied by the threshold behaviour of eq. (|3.2j ). This leads to ]T2 



f\ or \ 3 1 ft v m2 E n BW(s,F n ,M n ,T n ) 
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with 

rp2 A/f 5"p 

BW(s,F n ,M n ,T n ) = — " „ ■ (3.7) 

(M2 - s) 2 + M 2 r 2 n {s) 

M n , T n and F n are the mass, width and decay constant of the n-th resonance respectively. 
In the denominator we use an energy-dependent width T n (s): 

M n q(s 
r s q(M, 



r n (s) = — ^ r n , where q(s) = 7ry = yfts - s+)(s - s_) (3.i 



is the momentum in the centre-of-mass system. In the numerator, this energy dependence 
reproduces to the correct threshold behaviour of the spectral function and has been pulled 
out in writing eq. ( |3.6| ). Experimentally, the first two resonances, namely the _K"q(1430) 
and the _K"q(1950), are known and will be taken into account in our numerical analysis. 
All values for masses, widths and other parameters are given explicitly in section 5. 

The last quantity being required as an input for the phenomenological two-point func- 
tion is the scalar form factor at threshold, d(s + ). It can be calculated from an Omnes 
representation for d(s) || 



OC 



d(s + ) = d(0) exp j S -± f 5l{s,) ds' } , (3.9) 

71 J S'(S' — S + — IE) 

where Si(s) is the Ktt s-wave, 1 = 1/2 phase shift, which can be taken from experiment 
|35|-p7|. In writing eq. ( |3.9|) we have assumed purely elastic scattering. This is justified 



in the region close to threshold which dominates the integral. In section 5, we shall also 
compare with other methods to obtain d(s + ). 



3.2 The pseudoscalar two-point function 

In the pseudoscalar channel the lowest lying state with the relevant flavour quantum num- 
bers which contributes to the spectral function, is the i^-meson. In order to obtain good 
stability in the sum rule, we shall also include the next two resonances, for which there ex- 
ists experimental evidence: the K(14Q0) and the 7^(1830). For the i^-meson, a 5-resonance 
approximation is sufficient, because it is relatively long living, but for the two higher res- 
onances we again have to use an Breit-Wigner Ansatz with a finite width. As in the case 
of the scalar two-point function, we shall impose the correct threshold behaviour on the 
Breit-Wigner Ansatz. 

The next-higher intermediate state above the X-meson is the (i^7T7r)-system in an s- 
wave 1=1/2 state. Calculating the corresponding matrix elements for the divergence of 
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the axial-vector current to leading order in chiral perturbation theory, we obtain 
(n\d»(s W u)(0)\K + ( Pl )7r + (p 2 )n-(p 3 )) = 



(3.10) 



M 



K 



2 + 



(q + Pi)(p 2 - ps) + (q + p 2 )(pi - p 3 ) - (M| + m£ 



(n|^(s7 M 75«)(0)|K + (pi)7r (p 2 )7r (p 3 )) 



(3.11) 



M 



a' 



4 + 



(q + P2XP1 - Ps) + (<? + P 3 )(Pi " P2) - 2(M^ + M 2 



In our conventions F v = 132 MeV. 

Performing the phase-space integration, and modulating the threshold behaviour of 
the higher states with a sum of Breit-Wigner resonances, in the chiral limit the spectral 
function is found to be 

23 s En BW(s, F n , M n , T n ) ^ 



p5ph{S) 



F Z K M% S(s - M K ) + 



(3.12) 



3-2 12 tt 4 F* En BW(0, F n , M n , T n ) J 

The Breit-Wigner function is given in eq. ( |3.7|) . Notice that to leading order in the chiral 
expansion the threshold is shifted to s — 0. In this limit the width of the resonance becomes 
independent of the energy and, therefore, for consistency, we set T n (s) = T n . 



4 The perturbative continuum 



Above the energy sq, up to which experimental information on the spectral function is 
available, we approximate the remaining contribution by the perturbative continuum, ne- 
glecting all power corrections, being negligible for energies greater than sq. 
The two integrals which have to be calculated are 



I (u) = J e~ s/u p th (s)ds and h(u) = J ^e~ s/u p th (s) ds . 



(4.1) 



•50 



The expression for pth( s ) can De obtained from eq. ( [2.31) . Up to 0(a 3 ), we find 
N 



Pth{S) 



(M =F m) 2 s 1 1 + a (2c l2 L + en) + a 2 (3c 23 L 2 + 2c 22 L + c 2 i - c 23 7r 2 ) 



+ a- 



4c 34 L 3 + 3c 33 L 2 + (2c 32 - 4c 34 7r 2 )L + c 31 - c 33 vr 2 



(4.2) 



where L = ln(s//x ). Using the relations amongst the of eqs. (p,13j ) and ( B.14 ), one 
can convince oneself that also pth{s) satisfies a homogeneous RGE, as it should. 
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Hence, to calculate Io(u) and h(u), we need the following type of integrals: 

Q n—k 



I(a,n) = ( s"" 1 In" 4r e 
J a 2 

SO 



-s/u 



(da] 



n—k 



(4.3) 



with y = sq/u, and T(a, y) is the incomplete T-function [0^1 . Using this formula, we obtain: 

A' 



io(u 



8tt 2 



(M =F m) V T(2, y) 1 + ac u + a 2 (c 21 - c 23 ir 2 ) + a s (c 31 - c 33 n 2 ) (4.4) 



+ 2r'(2, y) ac 12 + a 2 c 22 + a 3 (c 32 - 2c 34 n 2 ) + 3r"(2, y) 



a 2 c 23 + a 3 c 33 



+ 4r" / (2,y)a 3 c 3 4 



The logarithms have again been summed up leading to running quantities evaluated at 
li 2 = u. The corresponding expression for I\ (u) can be obtained through replacing r( n ) (2, y) 
by T^{3,y) in eq. (O). 



5 Numerical Results 
5.1 Scalar two-point function 

In our numerical analysis of the sum rules, we shall mainly discuss the values of our input 
parameters, their errors, and the impact of those errors on the strange quark mass obtained. 
In the OPE, we generically have to evaluate expressions of the type a m (u)M n (u), where m 
and n are integer exponents. To achieve a systematic expansion in perturbation theory, we 
have decided to expand these terms consistently in inverse powers of L = ln(w/A|^). The 



relevant formula is given in eq. (|B.11| ) in appendix B. Besides the QCD scale parameter 
Ajfg, it involves the RG invariant quark mass M. However, as our analysis shows, the 
value for rh s depends strongly on A^r: it is roughly proportional to the leading term L. 
We therefore present our main results in terms of m s = m s (lGeV), which only displays 
a mild dependence on Aj^. In fig. 1, we show m s and rh s as calculated from the sum 
rule of eq. ( |1.9|) for A-^jg = 280, 380 and 480 MeV. The thick lines correspond to m s and 
the thin lines to rh s . Our Ajjg corresponds to 3 flavours and to NNLO, and has been 
chosen such that a s (Mz) = 0.118 ± 0.006, which covers most values obtained from recent 
analyses p9 |. All other parameters have been set to central values which will be discussed 



in the following. 

Choosing the stability region from which we determine m s to be in the range W = 
y/u = 2.0 - 3.0 GeV, from fig. 1 we obtain m s = 196, 189 and 198 MeV for A^ = 280, 
380 and 480 MeV. The lower end of the stability interval has been chosen such that the 
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O(a) correction is ^ 50% of the leading term, and at the upper end, the sum rule be- 
comes insensitive to the hadronic resonance structure and is completely determined by the 
perturbative continuum. The parameter so which determines the onset of the perturbative 
continuum has been adjusted so as to obtain optimal stability. For the three values of 
we find so = 5.5, 6.4 and 7.5 GeV 2 respectively. It is gratifying that y/so = 2.5 ± 0.2 GeV 
is just found to be around the energy at which the next resonance, which has not been in- 
cluded, is expected. However, because of the imperfection of our phenomenological model, 
we still have a residual dependence on sq. This dependence is displayed in fig. 2, where 
we plot m s for Ajjg = 380 MeV and s = 5.9, 6.4 and 6.9 GeV 2 . We do not include this 
variation in the error on m s because it is implicitly covered through varying all other input 
parameters to be discussed below. 

Let us next discuss the determination of the parameters in the phenomenological two- 
point function which has been presented in sect. 3.1. We have calculated these parameters 
from the experimental data for i^-scattering given by Estabrooks et al. [p?[] and Aston et 



al. |35|,|36[]. Generally, the s-wave amplitude and phase shift as and 4>s can be decomposed 
as follows: 

as = ~r ( Vie 2lSl - 1 ) + \. ( r/ 3 e 2 ^ ~ 1 ) , (5.1) 

where 771,3 and #1,3 are the elasticities and phase shifts for the I = 1/2 and 1 = 3/2 channel 
respectively. In writing eq. (|5.1| ) we have adopted the normalisation of refs. [p5} -|37[| for as- 
In ref. [fTTH , Si has been measured for yfs = 0.73 - 1.30 GeV, and S3 for 0.73 - 1.72 GeV. In 
this work it was also demonstrated that below 1.3 GeV, the scattering is purely elastic, that 
is, 771 = r]3 = 1. On the other hand in ref. |36| only as and </>s up to \/s = 2.52 GeV were 
given. In order to be able to calculate Si from the data by Aston et al. above 1.3 GeV, we 
have to subtract the 1 = 3/2 contribution. We fitted this contribution to a pure effective 
range 

tan5 3 (s) = aq(s)[l + f3q(s) 2 ] , (5.2) 

where q(s) is the centre of mass momentum defined in eq. ( |3.8|) . For the fit we used the 
full data set of ref. [37], however multiplying their error by a factor of 2 to get a x 2 /d.o.f. 
of order 1, and, below 1.3 GeV, for the data of ref. |36[ we calculated £3 from eq. ( |5.1| ) with 
Vi = V3 = 1- O ur best fit is obtained for: a = — 1.04GeV _1 and (3 = — 0.67GeV~ 2 . The 
value for a corresponds to a scattering length of —0.20 fm. 

Using this fit, we have then calculated 5\ for the data of both groups below 1.7 GeV. 
The result obtained for Si can be fitted to the sum of an effective range and a Breit-Wigner 
resonance: 

Si(s) = S E r(s) + Sbw(s) , (5.3) 
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where 5er is again of the form given in (|5.2|), and 



, c / \ M R T R (s) 

tan<WO) = —2 r. (5.4) 

{IV1 R — S) 



T R (s) has been defined in eq. ( |3.8| ). The lowest lying resonance in this channel is the 
#5(1430), and for our best fit we obtain: M K * (im = 1423 ± lOMeV, T K * {U30) = 268 ± 
25 MeV, a = 2.06 GeV" 1 and (3 = -1.37 GeV -2 . In this case a corresponds to a scattering 
length of 0.41 fm. Mk*(iaw) and IV*( 1430 ) are in good agreement to the values given 
in [35] . We show our fit in fig. 3 together with the data points and a fit to a pure effective 
range below 1.3 GeV. We do not give errors for a and (3 explicitly, because they are not 
direct input parameters, but we have varied them for the calculation of d(s + ). We have 



also compared our fit to 8% with other approaches, namely the /c-matrix formalism [40 



chiral perturbation theory [41, 42 and the quark Born diagram formalism [43|, finding good 



agreement. However, the x f° r our effective range plus Breit-Wigner Ansatz is lowest. 

Above 1.7 GeV, we do not know how to subtract the 1 = 3/2 contribution. We therefore 
assume two extreme cases for our evaluation of the errors on d(s + ). As a lower bound for 
Si, we use the pure effective range also displayed in fig. 3 up to infinity, and as an upper 
bound, above 1.7 GeV we take 5± = 180° constant ||. Using eq. (|3.9|) with 

d(0) = (M 2 K - Ml) f + (0) = 0.22 GeV 2 , (5.5) 

and varying a and (3 within the 1 a level, we find 

d(s + ) = 0.33 ±0.02 GeV 2 . (5.6) 

This value can be compared with other methods to calculate d(s + ). Close to s = 0, d(s) 
can be approximated by a linear function, 

d(s) = d( )(l + Ao^). (5.7) 



Using Aq = 0.019 f4"4], j4"5| , we get d(s + ) = 0.31 GeV . However, we know that d(s) is a 



concave function of s [45]. Thus the latter value should lie on the low side of the true value 
for d(s + ). On the other hand, calculating d(s + ) from chiral perturbation theory f45[ , we 
obtain d(s+) = 0.35 GeV 2 , also within the errors of our determination from experiment. 
Nevertheless, improved data for low energy #7r-scattering would be of great help for a 
more precise determination of the hadronic parameters. 

For the mass and the width of the second resonance #o(i950) we use the values found 
by Aston et al: M^ (1950) = 1945 ± 22 MeV and r^. (1950) = 201 ± 86 MeV. In addition, 
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we need the ratio of the decay constants for the first and second resonance which deter- 
mines their relative strength. From the dual-model vertex function, one finds the following 



parametrisation |TU 



£ - "4- 

In our analysis we shall use 7 = 0.5 ± 0.3. The lower value means almost no contribution 
of the second resonance, whereas for the upper value, the second resonance is nearly as 
strong as the first one. Both being extreme cases, we think that the above range for 7 is 
a conservative one. The induced uncertainty on m s from this parameter then turns out to 
be ± 6 MeV. 

Taken together, the errors in all phenomenological parameters lead to an uncertainty 
in the strange quark mass of ± 17 MeV, where we have added the errors quadratically. The 
dominant uncertainties in this set of parameters are due to d(s+) and rV*(i43o); namely 
± 12 MeV and ± 9 MeV respectively. A detailed table containing the relevant parameters 
entering the sum rule, their errors, and the resulting uncertainty in m s can be found in 
appendix D. 

For the theoretical part of the sum rule, we still have to discuss explicit quark mass 
corrections, the non-perturbative condensate parameters and the O (a 3 ) correction. As it 
turns out, because we work at a relatively high scale (2.0 - 3.0GeV), the quark mass 
and condensate contributions — except, of course, for the global quark mass factor - 
to the scalar sum rule are negligible. For completeness we nevertheless present our input 
parameters. For the light quark masses, we adopted the Particle Data values, m u = 
m u (l GeV) = 5 MeV and = m<f(l GeV) = 9 MeV For the quark condensates we 
assume SU(3) flavour symmetry, because anyhow their contribution is tiny, and take a 
standard value 

(qq)(lGeV) = - (225 ± 25 MeV) 3 , (5.9) 

where q = u, d, s. Let us, nevertheless, remark on the flavour dependence of the quark 
condensate. Because we use non- normal-ordered condensates, the strange condensate (ss) 
is closer to the light quark condensate (qq) than had we used normal-ordered condensates. 
Using eq. ( |C.1|) , we find ('.ss:)/ (ss) ~ 0.88. Together with the known value of the flavour 
breaking for the normal-ordered condensate ||[|7],[|8[ this leads to (ss) ~ 0.7 (qq). 



For the gluon condensate we use an average over most recent determinations, however 
with a quite conservative error, 

(aFF) = 0.024 ±0.012 GeV 4 . (5.10) 
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This value corresponds to 1 - 3 times the standard SVZ value |1|]. Finally, for the mixed 
condensate we take [fjgfl 



(gqaFq) = M 2 (qq) with M 2 = 0.8 ± 0.2 GeV 2 . (5.11) 

Again for the mixed condensate we can study its flavour dependence. Using the RG 
invariant combination of dimension 5 [p2| , |3~4| , we find (:gsaFs-) / (gsaFs) ~ 0.99. Therefore, 
the flavour breaking for the normal- and non-normal-ordered mixed condensates are almost 



equal. Using the values of refs. ]49],^(J, we thus have (gsaFs) ~ 0.75(gqaFq). 

The last uncertainty which we have to discuss, and unfortunately the dominant one, 
stems from the 0(a 3 ) correction. At W = 2.5 GeV and = 380 MeV, the 0(a) correc- 
tion amounts to ~ 40% and the 0(a 2 ) correction to ~ 30% of the leading term. Because 
the perturbative expansion has only asymptotic character, we conclude that the uncer- 
tainty coming from the as yet unknown 0(a 3 ) correction could be as large as 30%. In turn 
this would correspond to C31 ~ 580 if the correction is positive. For the error on m s we 
then find ± 27 MeV. This uncertainty turns out to be approximately halve of the error on 
the O (a 3 ) correction, because the sum rule scales like m 2 . Let us again emphasise that 
this is not completely consistent, because /3 4 and 74 are unknown too, and hence have not 
been included. From all of the above, our main result for m s from the scalar sum rule is: 

m s = m,(lGeV) = 189 ±32 MeV. (5.12) 



We shall, however, try to get some idea about the O (a 3 ) correction. From eq. ( |2.8| ) it is 
seen that the first 3 coefficients in the perturbative expansion grow almost geometrically. 
If this geometric growth continues for the 0(a 3 ) term, we would find C31 ~ 440. This 
would correspond to a negative 5% correction, leading to m s = 178 ± 18 MeV if the 
phenomenological error is included. 

So far, we only dealt with the sum rule of eq. (|1.9|) . Let us now briefly analyse the 1st 
moment sum rule ( |1.10| ). In fig. 4, we show the 1st moment sum rule for = 380 MeV 



and So = 5.4 GeV 2 in comparison to the 0th moment sum rule at so = 6.4 GeV 2 . All other 
parameters take their central values. It is obvious that the 1st moment sum rule is less 
stable. This is due to the fact that, because of the additional factor s in the dispersion 
integral, higher states are less suppressed and the sum rule becomes more sensitive to these 
corrections. Within the errors the value of m s thus obtained agrees with our previous result. 
Because the 1st moment sum rule is less stable, we refrain from taking an average with the 



value of eq. (5.12), but we take it as a corroboration of this result. 
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5.2 Pseudoscalar two-point function 

The hadronic input parameters for the higher resonances in the pseudoscalar channel are 
less well established than in the scalar channel. There are indications for two resonances 
with the same quantum numbers as the if-meson, the if (1460) and the if (1830) JKJ. The 
reported widths are rV( 1460 ) = rV( 1830 ) ~ 250 MeV, although the error is probably large. 
Since we shall show in the following that anyhow the pseudoscalar sum rule at present 
should be abandoned, we just use the central values. The ratio of the decay constants of 
the two resonances is again given by eq. ( |5.8| ) taking 7 = 0.5. 

For the input parameters on the theoretical side of the sum rule, we use the same set of 
values as for the scalar sum rule in the previous section. In fig. 5, we display m s calculated 
from the pseudoscalar sum rule for Ajj$ = 280, 380 and 480 MeV. The values for m s , 
obtained from this figure are 144 MeV, 136 MeV and 143 MeV for s = 4.9 GeV, 5.9 GeV 
and 6.9 GeV, and the three values of Ajjg, respectively. Even though the errors on the 
0(a 3 ) correction are large, we would expect that the central values for m s from scalar and 
pseudoscalar sum rules agree within the errors of the hadronic parametrisation, since the 
radiative corrections for flavour non-singlet scalar and pseudoscalar two-point functions are 
the same. This is, however, not the case. Let us comment on the origin of this discrepancy. 

Because we evaluate the sum rule at a rather high scale, 2-3 GeV, the dominant 
hadronic contribution comes from the if (1430). Therefore, the error on the global coeffi- 
cient of this resonance determines the largest error from the hadronic spectral function. In 
the case of the scalar sum rule, the corresponding quantity was the scalar form factor at 
threshold, d(s + ). Had we used instead the leading order chiral perturbation theory result, 
d(s + ) = d(0) = 0.22 GeV 2 , the central value of m s would have been 127 MeV — m s scales 
linearly with d{s + )\ The corrections from NLO chiral perturbation theory were large and 
led to d(s + ) = 0.35 GeV 2 , quite close to the value determined from experiment. 

For the pseudoscalar spectral function, we neither have available a NLO chiral per- 
turbation theory calculation of the (if7T7r)-contribution, nor is there enough data in this 
system to calculate the corresponding pseudoscalar form factor from experiment. Never- 
theless, an additional hint that the higher order corrections are large comes from a NLO 
chiral perturbation theory calculation of the matrix elements in if/4-decays f5l |, which are 
related to the matrix elements of eqs. ( ^.10[ ) and ( |3.11| ). We shall return to a calculation 
of the pseudoscalar spectral function at NLO chiral perturbation theory in a forthcoming 
publication. 

A different possible parametrisation of the pseudoscalar spectral function is the zero- 
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width approximation also for the higher resonances, 

p 5 {s) = F*M 4 K | 6{s - M 2 K ) + r K 6{s - M*) + r Kl 2 5{s - M'g) } , (5.13) 

where 

r * s K = 7±I < 5 ' 14 ' 

was also determined from a sum rule analysis ||. Here we adopted the abbreviations 
K' = K(1460) and K" = X(1830). If we insert this parametrisation of the spectral 
function into the sum rule of eq. ( |1.9| ), we find as a central value m s = 162 MeV, on the 
low side of the strange mass region as calculated from the scalar sum rule. Let us point out, 
however, that the zero-width approximation only serves as a lower bound on the spectral 
function, and the corrections are expected to be of the order of 10 - 20% [fTO] . In addition, 
it is unclear whether the error on tk given in eq. (|5.14 ) is a conservative one. 



To conclude, due to unknown uncertainties in the hadronic parametrisation of the 
pseudoscalar spectral function, at present, we shall not use the analysis as an independent 
determination of m s , but we plan to return to this issue in a future publication. 



6 Discussion 

Before we summarise our main results, let us discuss a few previous determinations of m s 
which are of interest in view of the present analysis. 

In 1982, Gasser and Leutwyler || reviewed the status of quark mass determinations 
known at that time. For the determination of m s , they used a combination of different 
QCD sum rules, together with quark mass ratios from chiral perturbation theory. Their 
final quoted value was m s = 180 ± 50 MeV. A calculation using finite energy sum rules 
at NNLO for scalar and pseudoscalar two-point functions was performed by Gorishny et 
al. ||. For the hadronic parametrisation of the two-point functions they assumed the zero- 
width approximation and found m s = 195 MeV, where the error was estimated to be of 
the order of 25%. 

The hadronic parametrisation of the two-point function in the scalar channel was im- 



proved in refs. H.^Q. 12] by taking resonances with a finite width. Whereas the authors 



of H only gave a lower bound on the strange quark mass, the final result of ref. JT2[ was 



m s = 192 ± 15 MeV. Since this was only a NLO calculation, due to the large radiative cor- 



rections, the error in this determination appears to be too small. In ref. |13j , Dominguez and 
de Rafael determined the strange mass by first calculating m u + rrid from the pseudoscalar 
sum rule for light quarks and then applying quark mass ratios from chiral perturbation 
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theory to find m s = 199 ± 33 MeV. They also refrained from using the pseudoscalar sum 
rule for the i^-channel directly because of large uncertainties in the hadronic parametrisa- 
tion. Although it is not obvious that the parametrisation of the spectral function is more 
reliable in the case of light quarks m u and ma, their central value is in good agreement 
to our results. The work of Narison is summarised in refs. H[l7j. By purely analysing 
the pseudoscalar channel, he obtained m s = 169 ± 10 MeV, whereas his final result for an 
average over determinations from different sum- rule channels gave m s = 160 ±9 MeV. The 
further reduction of m s is caused by the rather low value found in the (f) -meson channel, 
although the errors are very large. 

The last sum-rule determination which we would like to discuss, and the one which 



comes closest to our analysis of the scalar channel, is the work by Dominguez et al. [19|. The 
main differences are the following: they only used NLO expressions for the perturbative 
two-point function, they did not include the cancellation of mass singularities for the 
gluon condensate and the hadronic spectral function had been parametrised differently 
Nevertheless, their final result, m s = 194 ± 4 MeV, is in good agreement to our findings. 
However, the quoted error was obtained by varying solely Aj^s- Since this is a minor 
source of uncertainty, to our minds, their error is largely underestimated. Very recently, 
the strange quark mass was calculated from lattice QCD [[EJ. The result was fn s (2 GeV) = 
127 ± 18 MeV. Scaling this value down to \l = 1 GeV, we find m s = 180 ± 25 MeV. 

All determinations of the strange quark mass from QCD sum rules which have been 
discussed above are in general agreement with the calculation presented in this work, 
although in some analyses the errors were underestimated. Our main result was obtained 
from a sum rule for the divergence of the vector current: m s = 189 ± 32 MeV. In this 
determination, the error is dominated by the unknown perturbative 0(a 3 ) correction. We 
tried to estimate the 0(a 3 ) correction by observing a geometric growth of the first three 
terms in the perturbation series, and assuming that this geometric behaviour continues for 
the next order. The result then is: m s = 178 ± 18 MeV, where the error now dominantly 
stems from the hadronic parametrisation of the spectral function. Especially this latter 



value is in astonishing agreement with the lattice calculation of ref. ||52|| . The sum rule for 
the divergence of the axial-vector was found to be unreliable, due to large uncertainties in 
the hadronic parametrisation of the spectral function. Recently, instanton contributions 
to the determination of light quark masses were investigated |53| , and it was found that 
they can be sizeable if the scale at which the sum rule is evaluated, is too low. Because of 
the relatively large scale at which we obtain the strange mass, in the case at hand these 
contributions can be safely neglected. 

As a final application of our results, we estimate the masses of the light quarks m u and 
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ma from quark-mass ratios calculated in chiral perturbation theory, and the mass of the 



strange quark obtained in the present work. The relevant quark-mass ratios are [45,54 



' hl '" 25.7 ±2.6 and md m " = 0.28 ±0.03. (6.1) 



m u + m d m u + m d 

Using m s = 189 ± 32 MeV, we find m u = 5.3 ±1.5 MeV and m d = 9.4 ±1.5 MeV, whereas 
for m s = 178 ± 18 MeV we obtain m u = 5.0 ± 1.1 MeV as well as m d = 8.9 ±1.1 MeV, in 
good agreement to earlier determinations published in the literature. We plan to directly 
investigate quark-mass ratios from QCD sum rules in the future. 
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Appendices 



A The Borel transform 

The Borel operator B is defined by (s = Q 2 ) 



-s) n d n 



B u = lim 

s/n=u v / \ / 



The Borel transformation is an inverse Laplace transform pd^ . If we set 

oo 

f(u) = B u \f(s)}, then f(s) = f -f{u)e- s ' u du. 
L J J u 



We need the following Borel transforms JSC 

1 



-In" 4 



E( 

fc=0 



n— fc 



<9 



n—k 



{da) n ~ k T(a) 



and explicitly for n = 0, 1, 2, 3: 



1 



1 



By 



S a fi 2 



In 2 



1 i 3 S 

— hi 3 — 

s° /i 2 



u a T(a) 
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w a r(a) 
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u°T(a) 
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u a T(a) 



Here, L = lnw//i 2 and 



Z + V(a) 

Z 2 + 2^(a)L + ip 2 (a) - 

L 3 + 3^(a)L 2 + 3(?/> 2 (a) - ^(a))Z + </> 3 (a) 
-3^(a)^'(a) + <(a) 



For integer values of a, we have the following useful formulae 

n-1 



ijj(n) 



1 



i=i 



" n—1 



E^T-c(fc + i) 



i=i 



k > 1 . 
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To arrive at eq. (|1.6|), we also need: 



B„ 



(x + s) c 



u a T(a) 



-x/u 



(A.7) 



B Renormalization group functions 

For the definition of the renormalization group functions we follow the notation of Pascual 
and Tarrach |57fl, except that we defined the /3-function such that /3i is positive. The 
expansions of (3(a) and 7(a) take the form: 



f3(a) 
with 



0i 



Pia — (3 2 a 2 — j3 3 a 3 — . . . , and 7(a) = + 7 2 a 2 + 73a 3 + . . . , (B.t 



1 r 

6 



lliV-2/ 



Ik 



12 



17N 2 - 5Nf - 3C F f 



(B.9) 
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7i 



72 



48 



97iV + 9C F - 10/ 



(B.10) 



73 32 



U413 . iv2 _H£ iV C F -(^ + 24C(3))]V/ + ^Cj-(23-24C(3))C F /-|/ 2 



108 4 V 27 

where 73 has been taken from ref. [58| . 

We also need the expansion of an arbitrary product a m (u)M n (u) in inverse powers of 
L = ln(w/A 2 ): 



a m (u)M (u) 
with 



£j^n /2\ m + n 7i//3i 



(3 r { 1 \L 
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^10 = (A 72 - /?27l) > 

4/3 2 / \ 
d 21 = d u d w - —r± (mf3 2 + nj2) , (B.12) 
Pi 

1 Id 
- (d 21 - d n dx + d?„) + -3 (m/5 3 + n7 3 ) - (/?i/3 3 - Pi 



■ * -$ 2 A 

2 a n ^2 Gl1 ' 
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Here, M is the renormalization group invariant quark mass. Because 7i//3i is close to 1/2, 
from the prefactor of eq. (p.ll| ) we read off that two powers of M count approximately 
like a in the perturbative expansion. This way of counting orders in perturbation theory 
has been adopted in this work. 

Using the renormalization group equation (|2.5|) , we find the following relations between 
the c,- 



Cl2 

C22 
C23 



2 7lC01 



+ 27i)cn + 272C01 



— 71 (/3i + 27i)c i 



(B.13) 



C32 
C33 
C34 
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4 
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6 
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(& + 272)012 + 2(/3i + 7i)c 22 



(B.14) 



C Renormalization group invariant condensates 



To the order considered in this work the renormalization group invariant condensates are 
given by |31"|. |3~4|1 : 

1 



where 



+ w 2 



mi rrq 



fi4 = {aFF )^ 4 faJ2mi(qiqi)-fJ2mt, 
Pi i Pi i 



Wi 



7? 



w 2 



I(4 Tl -/5 1 )--^(4 T2 -/5 2 ) 
6 471 



(47i -A)' 

and 7° = N/(2tt 2 ) is the leading order anomalous dimension of the vacuum energy. 



(C.l) 
(C.2) 

(C.3) 
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D Input parameters and the error on m s 

We have omitted those input parameters from the table whose variation results in an error 
on m s less than 1 MeV. 



Parameter 


Value 


Am s 




380 ± 100 MeV 


±5 MeV 


0(a 3 ) 




±27 MeV 


d(s+) 


0.33 ± 0.02 GeV 2 


± 12 MeV 




1423 ± 10 MeV 


±3 MeV 


IV* (1430) 


268 ± 25 MeV 


±9 MeV 


rV*(1950) 


201 ±86 MeV 


±4 MeV 


7 


0.5 ±0.3 


±6 MeV 
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Figure 1: m s and m s for three different values of A^r. Thick lines correspond to m s and 



thin lines to fh s . 
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Figure 2: m s for three different values of s and = 380 MeV. 
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Figure 3: Our fit for the s-wave 1=1/2 phase shift, compared to the experimental data. 
The thin line corresponds to a pure effective range. 
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Figure 4: Comparison of m s from the Oth and lst-moment sum rules. 
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Figure 5: m s calculated from the pseudoscalar sum rule. 
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